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Ultra-high energy [E^, > 10® GeV) neutrino-nucleon total cross 
^ ■ sections a'^^ are estimated from a soft non-perturbative model com- 

• plemented by an approximate QCD evolution, explicitly calculated. 

The resulting asymptotic energy behavior of the neutrino-nucleon cha- 
rged-current total cross section is derived analytically and predictions 
concerning the observation of ultra-high energy neutrinos in future 
experiments are presented. 
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1 Introduction 



Calculations of ultra-high energy (UHE) neutrino-nucleon cross sections 
[Ill21IHlllllSlinilZllH] have attracted attention for many years since, on the one 
hand, UHE neutrino fluxes are predicted by cosmological models but have 
not yet been measured experimentally and, on the other hand, theoretical 
models exist for the evaluation of these cross sections. 

Neutrino telescopes are the most promising tools for probing the distant 
stars and galaxies. Due to their high stability and neutrality, neutrinos arrive 
at a detector on a direct line from their sources, undefiected by intervening 
magnetic fields, with expected energy up to 10^° eV. Whereas high-energy 
photons are completely absorbed by a few hundred grams/cm^ of material, 
the interaction length of a 1 TeV neutrino is about 250 kilotonnes/cm^ jS]. 

An important scientific goal of large-scale neutrino telescopes is the de- 
tection of UHE cosmic neutrinos produced outside the atmosphere, i.e. neu- 
trinos produced by galactic cosmic rays interacting with interstellar gas, and 
extra-galactic neutrinos. High neutrino energy brings a number of advan- 
tages [SI El- Firstly, the charged-current cross section increases as a ~ 
for Ey < 10^^ eV and is believed to rise as cr ~ for higher Ey. (In this 
paper we question the last point.) Secondly, the background of atmospheric 
neutrinos decreases as compared to the flux from extra-galactic sources, ap- 
proximately as E~^-^. 

Another important and promising aspect of UHE neutrino physics is the 
possibility to study the production of exotic objects, such as supersymmetric 
particles P or black holes 0. 

The phenomenology of UHE neutrinos and their detection depend on 
the neutrino-nucleon total cross section, which has been calculated in the 
standard model (see e.g. ^3j). A striking feature of these calculations is 
the continued power-law growth of the cross sections with E^ at highest 
energies. The rise with E',^, on the other hand, is directly related to the 
very low-x behavior of the nucleon parton distributions derived from e.g. the 
HERA data. Problems related to the unitarity bound have been discussed 
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in inmnnnms- 
At ultra-high neutrino energies the total cross section is completely dom- 
inated by its deep inelastic component. The energy dependence of neutrino- 
nucleon cross section, as suggested from the small-x nucleon structure func- 
tions (SF) measured at HERA, should show a rapid increase with energy, 
like pB], where A scales as ln((5^/A^) and ranges approximately in the 
interval 0.1 ^ 0.4 for 1 GeV^ < < 100 GeV^ (see Ref. HI). 

The highest available ep energy at HERA, y/s^ =314 GeV, is still too 
low as compared to the expected UHE neutrino-nucleon collision energies. 
For the neutrino energy ~ 10^2 Qgv, the relevant x- region is x ~ 10 , 
while the HERA measurements extend now to a; ~ 10~^; however such low 
values of x are measured at < 0.1 GeV^. For the interesting region of 
q2 ^ M^, SF are measured only up to a; ~ 10~^ in the DO experiment of 
inclusive jets [1^1. Since studies of the UHE uN cross sections began nearly a 
decade ago [H El El HEj our knowledge about quark distributions and nucleon 
SF improved significantly [T7] . 

In the present paper we study the UHE behavior of the neutrino-nucleon 
charged- current cross section starting from models where nucleon SF have 
a milder increase. Our analysis relies strongly on recent progress and our 
previous experience in the studies of the nucleon SF [T7 | IT ^ IT? H l ^ m i l ^ l^. 
based on phenomenological fits, extended by QCD evolution. 

The paper is organized as follows: in the next Section we briefly review the 
steps involved in the determination of the neutrino-nucleon cross section and 
illustrate our approach; in Section 3 we study the asymptotic behavior of the 
total neutrino-nucleon cross section; in Section 4 we discuss the implications 
of our results for the observation of UHE neutrinos in future experiments 
and finally in Section 5 we draw some conclusions. 
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2 Nucleon structure functions and neutrino- 



nucleon cross sections 

In order to compute the cross section for neutrino-nucleon interactions at 
high energies, one needs a detailed description of the quark structure of the 
nucleon. The special challenge of the UHE neutrino-nucleon cross section is 
that the VT-boson propagator emphasizes smaller and smaller values of x as 
the neutrino energy E,y increases. In the UHE domain, the most important 
contribution to the uN cross section comes from x ~ M^/{2MEi^) and 
hence the greatest uncertainty in the predictions comes from the small-x 
extrapolation. 

Deep inelastic scattering of UHE neutrinos is unique in the sense that 
it explores extreme regions of the phase space where no accelerators 

data exist. One thus must rely on predictions or assumptions regarding the 
nature of these distributions, in particular in the region of ultra-low x and 
low Q^, as well as in the region of low x and low Q^. 

Let us consider the double differential cross section for the deep inelastic 
scattering of neutrinos on an isoscalar target for charged-current processes. 
We will limit ourselves to the leading order corrections to the simple parton 
model and hence all parton model formulas remain unchanged except that 
the parton distributions now depend on x and and not only on x. In 
particular, the Callan-Gross relation, F2 = 2xFi or = 0, holds in leading 
order [21]. With the usual notation we can then write for charged-current 
neutrino interactions 



For anti-neutrino charged-current processes one must change the sign in front 
of F^{x,Q'^). (The changes necessary in order to describe neutral-current 
neutrino interactions are given in many textbooks and review articles, for 
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example in j25). In Eq. (Q), Gp is the Fermi constant, M is the nucleon 
mass and the variable y is related to the Bjorken x through the relation 

^ x{s -M^) xs ' ^ ^ 

where s is the square of the cm. energy for the neutrino-nucleon scatter- 
ing. The laboratory neutrino energy = {s — M^)/(2M) is approximately 
s/{2M) in the energy region of interest. 

The main purpose of this paper is to evaluate the asymptotic behavior of 
the total cross section for charged-current neutrino-nucleon process 



^0 



dx dy 



in a model for the structure function F2 based on a soft non-perturbative 
input complemented by an approximate QCD evolution. The contribution 
of xF^ can in fact be neglected at ultra-high energies. Unfortunately, exper- 
imental data for the neutrino SF do not extend to the "deep sea" region, 
i.e. to X smaller than 0.001, and to large Q^. Both these limits, for the 
electromagnetic structure functions F^p, are subject of debate regarding the 
presence of possible saturation effects at small x and discrepancies in the 
evolution for large Q^. F2 differs from Fj^^ at small as can be easily 
seen by formally integrating the double differential cross section (0) over x 
at fixed value of the product {E^y) 

1 da GlM 



lim — — = ^— / dxF^{x,Q' = 0) , (4) 
y^o E^ dy n ' 



formula that has been used for the normalization of the neutrino flux 
The integral in Eq. when evaluated for F2^, would vanish by gauge invari- 
ance. At increasing Q^, this difference will disappear but, at the moment, no 
smooth transition has been found in all the range of the variables appearing 
in Eq. ©. 
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2.1 The non-perturbative input 

In a number of papers [11111312011^112212312^1 a diffraction model lias 
been developed, describing equally well hadronic reactions at high energies 
and DIS at small x. The main idea behind the model is that SF at small 
and moderate values of are Regge-behaved. The Q^-dependence was 
introduced in a phenomenological way. 

It was shown [THl dl 1211 123 1211 ^^^^ data on deep inelastic ep scattering 
from HERA and elsewhere can be fitted perfectly well by a "soft" Pomeron 
alone. This Pomeron could be a dipole with a ln(l/x) dependence or a 
squared logarithm, saturating asymptotically the Froissart bound. For the 
electromagnetic SF, parametrizations exist ^71 that allow for a simplified 
evolution until very large Q^, while maintaining a soft dependence on the 
Bjorken x. As an example we will consider a structure function F2{x, Q^) of 
the form 

F!^{x,Q') = 

and derive for it an explicit approximate evolution starting from = Q^. 
The main motivation for our choice of the model is that at small there 
is similarity between the small-x behavior of the structure function and the 
s-dependence, x ~ Q'^/s, of the hadronic total cross section at high energies. 
The alternative approach, with a power behaved dependence, has been dis- 
cussed in the literature in several publications. The differences between the 
two approaches can be significant and will be discussed in Section 3. 

We notice that the coefficients a(0), 6(0) and i/(0) are constrained by 
Eq. (01) whose left hand side can be determined from the experiment, Ref. 
The constraint is 

where 7 denotes the Euler's constant and ip[z] is the logarithmic derivative 
of the F-function. When low-Q^ data will extend to lower x values, it 
will be possible to check the choice ©. 



a{Q') + b{Q')\n[- 

X 



X 



(5) 
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By imposing a lower cut in Q^, Q'^ = Ql, and rewriting Eq. Q as 

2ME^Jq2 Jq2/s X \dxdy 

where in the differential cross section y = Q'^/{xs), the F2 contribution to 
the total cross section can be written in the form 

uN I PN 
-uN _ " ^ " 



X 



2.2 evolution 

Let us consider a simplified approach to the DGLAP |2H| evolution, where 
the parton distribution functions satisfy the relation q[x) = q{x) = u{x) = 
d{x) = s{x) = 2c{x) = 2b{x) and F2, for five fiavors, is given by F2 = 8xq{x). 
This assumption gives Fg^^ = 17xq{x)/9 and, as shown in with this 
simple form for the SF a high quality fit to the HERA data can be obtained. 

For each parton distribution we assume the same function (jSj) chosen for 
F^{x,Q') 



ai{Q^) + bi{Q^)ln{^ 



(l-xr(«'), it = q,g) (7) 



with the same functional form for quark and gluon but with different coef- 
ficients. Experimental data for quarks and global analysis, for example the 
parton analysis of MRST J7] , for the gluon could determine the coefficients 
in the expressions for xfi{x,Q^), where i = q,g stands for quark or gluon, 
respectively. 

In the leading order of perturbation theory, an approximate evolution is 
feasible by using the results of Ref . [22] • The only condition is that ^ bi 
{i = q,g), where = ai{Ql) and bi = h^Ql), so that no interference appears 
in the evolution of the coefficients multiplying the different powers of the 
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logarithm^. We define /3o = H — 2//3, where / = 5 is the number of flavors 
in our approach, 



t = In 



am 



In 



ln(QVA^ 



and recall the form of the scaling variables, 



cr = 2\ —dggt ln(l/x 



p 



-d f 



a 



ln(l/x) 21n(l/x) 



(8) 



where dgg = —12/j3o ^ . Then, from [21], we get for small x the evolution of 
the structure function F2(x, Q^) = 8xfq{x,Q'^) = 8x[/~(x, Q^) + fq{x,Q^)] 
as 



xfg (x, Q 



ttg + bg 



(a+p/i((T) + 6+/o(a)) -e- 



+0(p), (10) 



where = 16//(27/3o), = 1 + 20//(27/3o), a+ = /(a^ + 4aj9)/9 

and an analogous expression for 6^. Ci(z^) is related to the logarithmic 
derivative of the F-function: Ci(z/) = ^/'[l + z/] — Io{z) and h{z) are 

modified Bessel functions of the first kind. 

We notice that the presence of higher powers of the logarithm of 1/x can 
be easily handled in this formalism provided that Oj ^ bi ^ .... While 
the evolution given by Eqs. (P)- ()10|) is correct in the limit of very large a 
and ^ Ql, the requirement that the initial input must be reproduced for 

= Qq introduces a corrective term in that does not affect the evolution 
for large Q^. Eq. (fTUjl should then be replaced by 



xf^ix, Q') = {a^phia) + 6+ [/o(a) + A]) 
where A = bgCi{u)/b+ - 1. 



-d+{i)t 



0{p) 



^The requirement ^ hi is needed for all values. It drastically simplifies our 
formulas since it avoids a mixing between and hi parameters during the evolution. 

^We follow here the notation of d is the ratio between the anomalous dimension 
and twice the QCD /3- function in leading order, while d and d are, respectively, the singular 
and regular part of d when n — > 1. 
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3 Asymptotic behavior of a 



uN 



We rewrite Eq. (jH)) in the form 



dx 

X 



-UN 



1+1 



2Gl 



xs 



dQ 



{xf-{x,Q') + xf, 



[X 



(12) 



where xf^{x,Q'^) are given in Eqs. (j9|)- |T0|) respectively. The integral over 
X can be done exactly for /'{xjQ"^). Let us denote this integral with the 
symbol J_ (s, Q^) . The result is presented in Eq. ()A.1|) of the Appendix, and 
the leading behavior, for large s, of the x-integral is given by ln^(s/Q^). 

The presence of many scales (A^, Ql, and s) makes difficult any pre- 
cise estimate of the asymptotic behavior of the xf'i^x^Q"^) contribution to 
a'^^ . From the Ostrowski's inequality [SO] for integrals it follows, however, 
that this contribution has an upper bound for s — >^ cxd: 

2 

12 



dQ' 



+ Mf 



w 



J-{s,Q') 



< Kin' — T 



(13) 



where K is a suitable constant. 

Only the small-x evolution of fq{x, Q"^) is known from pO|) and the con- 
dition p <^ 1 requires a cutoff at the upper limit of the x-integral. Let us call 
it xm with xm < 1, while the Q^-integral has now xms as the upper limit 
of integration. The first integral we meet, when considering the term with 

/+(x,g2),is 



)2/s ^ 



-d+{l)t 



(14) 



and, with the change of variable x = exp[—z'^/{—Adggt)], the solution can be 

easily found. Setting (Ju = 2\J ~dggthi{s / Q"^) and an = 2\J —dggt ln(l/a;Af), 
we obtain: 

Xi(s,Q2) = / dzh{z) = /o((T„) - /o(a,) (15) 
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and 



1 Z"'^" 1 / \ 

l2{s,Q^) = — ^ — / dzzIo{z) = — ^ — {auhicru) - (rili{(re)j . (16) 

Other integrals, coming from powers of Q^/{xs) in Eq. (fT^ . are 0{p) with 
respect to those in Eqs. (fT3j) and (fT^ and it is consistent with our approxi- 
mation in Eq. (jlOp to neglect them. The proof can be found in the Appendix. 
We define now the integrals 

Ms) = j^^ dQ' [ q^/m^ ) e-'-'^'^%is, Q'), {k = 1, 2) (17) 

and notice that their large-s behavior determines completely the contri- 
bution to the mean total cross section in this limit. Of the two integrals 
J^iis), 1/2(3) only the first one needs to be evaluated explicitly since 

Ji[s) + O 



dins \\ns 

We consider then the integral J7i(s). The function under the integral sign in 
J7i(s) vanishes at both integration limits and has a bell-shaped form with a 
well defined maximum. The maximum moves at larger values of Q^, when s 
increases, and its position is defined, in the limit s 00, by the solution of 
the transcendental equation 



(3oasiQl)Hs/Q'' 



lt\n{s/Q^) 



X ( IT^^e* + Poas{Ql)d+il) ) - 47re* = , (18) 



where the variable t has been rewritten as 



t = In 



1 + (^^^^ HQVQl) 



Eq. (fTHj) can be solved for ln(s/(5^), obtaining an explicit expression for s 
as a complicated function of the value of at the maximum Q^ax- ^® 
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rewrite Eq. (fTH|) as /(s, Q"^) = 0, we get the slope 



dQl 



max 



ds 



ds 



max 



that, as can be easily shown, vanishes when s ^ oo. 

In order to estimate the position of the maximum as a function of s, 
in the range 10^ GeV < s < 10^^ GeV, we have chosen the parameters in 
intervals that include their most probable value and found that the position 
of the maximum can be roughly reproduced by a term linear in (Ins)^'^, 
while the value of the integrand, at the maximum, is near the value of the 
Bessel function there. The width at half height of the peak goes like In s. 
By multiplying the value, at the maximum, of the integrand in Ji{s) by the 
width of the peak, we get the right order of magnitude for the corresponding 
integral in Eq. ()17p. as can be verified by a numerical computation. In this 
limit the mass of the W-hoson plays little role since Q^ax ^ until 
s 10^^ GeV^. It is then possible to conclude that, in this approach, the 
total cross section for neutrino-nucleon scattering increases with s faster than 
any power of In s but slower than any power of s. 

Most of the existing models of F2 assume a fast, power- like increase in 
l/x that transforms in a like increase of a''^ {E) [HElinillliniElIZllinilSl]. 
In an alternative approach [IHl HI 1201 IHl 1221 121 El the HERA data are 
shown to fit a model of the SF, at low and intermediate with a moderate, 
logarithmic, increase in l/x, the observed "HERA effect" being attributed to 
the decrease of the non-leading contributions, relevant at larger x. DGLAP 
evolution transforms the logarithmic behavior in the way we have shown 
in Eqs. (P)- (fTTH) . It is important, however, to realize that this asymptotic 
behavior is quite different from the one obtained starting with a power, x~^. 
According to fnil2niEHl! if the starting distributions have a power-law form 
in x, the leading term keeps the same x-dependence under DGLAP evolution 
also at the next-to- leading order. The component xf^{x, Q^) in Eq. (fTn|l can 
be represented, in the leading order of perturbation theory and in the small-x 
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limit, by the approximated form [^21 El] 

x/+(x,Q2) ^ Z^x/,(x,Q^)e-'^+* , (19) 

if xfi{x,Ql) = K{Ql)x-^, {i = q,g). In Eq. (HHD, d+ = 7+/(2/3o) is the 
largest eigenvalue of the anomalous dimension matrix evaluated at (1 + ^ 
For /i ~ 0.4, as in Ref. [THI, is negative and ^ —3.4. In this approach 
it is easy to show that a"^ oc {s/M^Y + 0(s°). 

It is interesting to compare our result based on Eqs. (fTHjl and (fTI)|) with a 
model where the cross section increases like a power of s . We notice that 
Eq. (fT^ can be written as 

^ ^ {a^Ji + b^J2) (20) 

and that all the integrals can be evaluated numerically, choosing for example 
Ql = 10 GeV^, asiQl) = 0-23 and the values xm = 0.001,0.01. The choice 
of Ql = 10 GeV^ is reasonable because in the region of near this value 
the rise in x of SF is well described by a ln(l/x) and/or ln^(l/x) behavior ^. 

By evaluating the integrals numerically, a comparison can be done with 
the values of a'^'^ obtained from Table I of Ref. j3] . We plot in Fig. ^the values 
of Pi = {4:G'jp/n)J'i and the result of Ref. [3]. As can be seen from Fig. [H the 
power-behaved cross section increases with s faster than the components of 
a'^^ in the relation pO|l . The explicit values of the parameters and 6^ do 
not enter this comparison, but it turns out that their order of magnitude can 
be made consistent with the parametrizations of the HERA data studied in 
Ref. [H]. 

■'Note that the power-hke behavior ~ x^'''"* of parton distributions is in agreement also 
with the DGLAP dynamics in the range near the starting point Qq (see Refs. |35II34| '). 
However, the power-hke behavior cannot be combined together with (|19|) to some unified 
Regge-hke Q^-evolution of parton distributions. 

''For larger values, the rise of F2 is stronger and is better described by a power 
law x-^, with A - 0.3 at - 100 GeV^. Evidently, our initial conditions ((TJ cannot 
be applied here. However, after evolution in our approach, we get a behavior with a 
Q^-dependent effective slope X^q{Q^) in agreement with experimental data |14| . 
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The Bessel-like asymptotic behavior of the total cross section in Eq. (|2Uj) 
violates the Froissart bound since it is steeper than any power of In s, though 
being flatter than any power of s. The comparison with a power-behaved 
cross section, strongly violating the Froissart bound, makes us confident on 
the fact that additional contributions needed to restore unitarity should be 
less relevant in our approach. 

4 Observation of ultra-high energy neutrinos 

The evolution in leading order considered in the present paper overesti- 
mates the cross sections. Next-to-leading calculations show (see, for exam- 
ple, Ref. [HI) that the ratio ctnlo/o'lo decreases when the neutrino energy 
increases and it is reasonable to assume that our estimates constitute an up- 
per bound for the true predictions of the model. Hence the NLO evolution of 
the soft non-perturbative input could lead to a more pronounced difference in 
the predictions with respect to a power-like input. Another possible cause for 
the decrease of the true cross sections, that has not been taken into account 
here, resides in the gluon recombination at small x. ^ This recombination can 
be mimicked by considering a large, "effective" gluon anomalous dimension, 
with a consequent increase of the value of c/+(l) that appears in the expo- 
nential exp(— t) in Eq. (fTUj). Despite these possible corrections, it is of 
interest to consider the constraints imposed by our model on the interaction 
length of neutrinos and compare them with a more conventional approach. 

For this purpose, we considered the (water equivalent) interaction length, 
defined as 



where A^^ = 6.022 x 10^^ cm ^ (water equivalent) is the Avogadro's number, 
and evaluated the ratio of Lint(-E'i/) to Lint(10'^ GeV) as obtained from our 

^The results of numerous analyses of the higher- twist contributions in DIS [HEI and of 
recent papers discussing screening effects in neutrino- nucleon cross-section |lllll2j are still 
not conclusive. 




a' 
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model, with /a^ ranging between 1/20 and 1/3 and setting Ql = 10 GeV^, 
<^s{Qo) = 0.23 and xm = 0.01. In Fig.Elour result is compared with the same 
ratio determined from the results of Ref. j^. At higher energies Lint in our 
model could be larger than expected in more conventional approaches. 

Experiments are planned to detect UHE neutrinos by observation of the 
nearly horizontal air showers (HAS) in Earth atmosphere resulting from u- 
air interactions. The expected rates are proportional to the neutrino-nucleon 
cross section. The resulting cross section at lO^'' eV is rather high, ~ 10~^^ 
cm^ according to the estimates of jlHl 0] , and results from the extrapolation 
of parton distribution functions far beyond the reach of experimental data. 

If the cross section is lower, then the event rate for neutrino-induced HAS 
is reduced by the same factor. This reduction would compromise the main 
detection signal that has been proposed for UHE neutrino experiments. A 
smaller cross section would however offer a double advantage for the planned 
experiments jl]: 

1. with a new search strategy, based on the observation of upgoing air 
showers (UAS) initiated by muon and tau leptons produced by neutri- 
nos interacting just below the surface of Earth, the neutrino event rate 
with a small cross section is actually larger than the HAS rate with a 
large cross section; 

2. the future detectors can also measure the neutrino cross section at ener- 
gies far beyond those achievable in collider experiments thus providing 
important information for particle physics. 

If the proposed model is correct, the above scenario may have better 
chances to be realized. 

5 Conclusions 

UHE neutrinos became an extremely rich and interesting testing field for 
various physical phenomena, such as cosmology, astrophysics (origin of UHE 
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neutrinos), theory of gravitation (production of black holes), unified theo- 
ries, testing the standard theory and beyond, as well as the strong inter- 
action theory, responsible for the hadron structure. We predict that the 
UHE neutrino cross sections rise moderately with energy, in accord with the 
above-mentioned "soft" Pomeron dynamics. 

In the present paper we concentrated on the calculation of the asymptotic 
behavior of the structure functions and resulting cross sections. We are 
aware, however, about the important role of the large and intermediate x 
behavior of the structure functions to be fitted to the existing, relatively 
low-energy, data. The cross sections calculated in this way should match 
with their UHE asymptotic behavior. We intend to address this point in a 
forthcoming publication. 
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[GeV] 

Figure 1: The curves of Pi (solid for xm = 0.001, dashed xm = 0.01) and P2/IO 
(dotted for xm = 0.001, dot-dashed for xm = 0.01) are plotted versus the neutrino 
energy E^, together with data for a'^'^ (diamonds) from Ref. 0. 
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[GeV] 

Figure 2: Li^tiEu) / Li^t{'^^^ GeV) versus in our model (solid lines) and ac- 
cording to the results of Ref. ^ (diamonds). 
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A Appendix 



In this Appendix we include formulas that are not strictly necessary to un- 
derstand the logical lines in the main text, but can be useful in reproducing 
some relevant results. 

In order to determine the asymptotic behavior of a'^^ we needed some 
cumbersome integrals and the first one is 



1 



^ dx 

Q^/s X 



1+1 



xs 



2\ 2' 



Xfq {X, 



5a„ 



3Ci(z/) 



40^ 



S2 



+ 21nM^)6, 



-S+ Ci(u) + - k 



(A.l) 



Next, we prove that the integral 



Q' r' dx ^. . 



2/„ X^ 



(A.2) 



QV 



is non- leading, 0{p), with respect to the integral evaluated in Eq. (fT3j) . The 
identity 

^ 1 

(A.3) 



^ [p"/„(a)] = -lp«+i/„+i(a), n = 0, 1,2, 



can be easily obtained by using the definitions (jH)). Integrating repeatedly 
by part, we get 



J X'^ X J X"^ 

^ oo 



(A.4) 



k=l 



The dominant contribution to the integral ()A.2|) comes from the lower limit 
of integration and, since the z asymptotics of Ik{z) does not depend on k, 
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we get 




s 



that is non- leading with respect to the integral in Eq. (|T5|) . By using the 
same arguments, it is easy to show that 
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